function y = Rs(x,k)
    y = k*(exp(-k*abs(x)))/2;
end

% 选择k值和x值的范围  
k = 1; % k=1  
x = linspace(-5, 5, 1000); % 生成从-5到5的1000个点  
% 计算y值  
y = Rs(x, k);
% 绘制图形  
figure; % 创建一个新图形窗口  
plot(x, y); % 绘制曲线  
title('Rs(x, k) with k = 1'); % 添加标题  
xlabel('x'); % x轴标签  
ylabel('Rs(x, k)'); % y轴标签  
grid on; % 显示网格

% 参数设置  
T = 10; % 时间范围（或空间范围），足够大以展示函数的全貌  
dx = 0.01; % 采样步长  
x = -T/2:dx:T/2; % x向量  
  
% 生成信号  
signal = Rs(x, k);  
  
% 计算FFT  
N = length(signal); % 信号长度  
X = fft(signal); % FFT结果  
  
% 计算幅度谱  
magnitude = abs(X); % 幅度  
magnitude_normalized = magnitude / N; % 归一化幅度 
  
% 频率向量（对于FFT结果）  
fs = 1 / dx; % 采样频率（假设dx是采样间隔）  
f = (0:N-1) * (fs / N); % 频率向量  
  
% 注意：由于信号是实数的，FFT结果是共轭对称的  
% 我们通常只绘制前N/2+1个点（不包括0频率的直流分量，如果你想要包括它也可以）  
  
% 绘制幅度谱  
figure;  
plot(f(1:round(N/2)+1), 20*log10(magnitude_normalized(1:round(N/2)+1))); % 转换为dB并绘制  
title('Magnitude Spectrum of Rs(x, k)');  
xlabel('Frequency (Hz)');  
ylabel('Magnitude (dB)');  
grid on;  
xlim([0 fs/2]); % 只显示正频率部分  
  
% 注意：这里的频率单位是Hz，但实际上是基于采样率和FFT索引的。  
% 如果dx的单位不是秒，则fs的单位也不是Hz，但概念仍然相同。
